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Abstract 

We study the magnetization dynamics of a molecular magnet driven by static and variable 
magnetic fields within a semiclassical treatment. The underling analyzes is valid in a regime, 
when the energy is definitely lower than the anisotropy barrier, but still a substantial number 
of states are excited. We find the phase space to contain a separatrix line. Solutions far from 
it are oscillatory whereas the separatrix solution is of a soliton type. States near the separatrix 
are extremely sensitive to small perturbations, a fact which we utilize for dynamically induced 
magnetization switching. 

PACS numbers: 
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I. INTRODUCTION 



Molecular magnets (MM) are molecular structures with a large effective spin (S), e.g. 
for the prototypical MM Mn\2 acetates [1] S — 10. MM show a number of interesting 
henomena that have been in the focus of theoretical and experimental research [l, 2, 3, 4, 



q y, y, q, 



To name but few, as a result of the strong uniaxial anisotropy, MM 
show a bistable behavior jlj]; they also exhibit a resonant tunnelling of magnetization 
that shows up as steps in the magnetic hysteresis loops 

aaa. 

Of special relevance for 



applications in quantum computing is the large relaxation time of MM 

The present theoretical work focuses on the dynamics of the magnetization. The 
established picture of macroscopic quantum tunnelling of the magnetization is as fol- 
lows: The MM effective spin Hamiltonian H = — DS 2 possesses degenerate energy levels 
±Mg, —S < M s < S separated by the finite barrier Eb = DS 2 . At low temperatures only 
the lowest levels Ms = ±5* are populated. Those two states are orthogonal to each other and 
no tunnelling is possible. An anisotropic perturbation E(S 2 — S 2 ) does not commute with 
the Hamiltonian H = —DS 2 and mixes therefore the states at both sides of the anisotropy 

n 

barrier leading thus to tunnelling |8j. Reversal of the magnetization due to macroscopic 
quantum tunnelling has a maximum for the states close to the top of the barrier. This case 
corresponds to the high temperature limit. 

In this work we consider the magnetization dynamics induced by constant and harmonic 
external magnetic fields: The influence of a variable magnetic field on MM at low tem- 
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111 ] . It was shown 



peratures, i.e. when only 2-3 levels are excited was considered in 
that in this case the problem is reduced to a three level Jaynes-Cummings model, the so 
called Lambda configuration. Therefore, it is analytically solvable in principle. The low 
temperature assumption is however quite restrictive 12j, |l3j: If only the levels E°, E 1 , E 2 



are involved the low temperature approximation is applicable for temperatures obeying [11] 
ksT < E 1 — E°, where ks is the Boltzmann constant. For Mn\2 this leads to the esti- 
mate T < Q.QK [3]. Obviously, if the temperature exceeds T, an approximation with a 
large number of levels participating in the process is more appropriate. In this case the 
quasi-classical approximation for the spin dynamics becomes applicable 
aim here to conduct such a study. MM will be modelled as in previous studies, e.g. in 
Refs. 10, [ll[ . We consider the dynamical reversal of the magnetization, caused not by an 



Q Q. It is our 
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anisotropic perturbation but by a constant and varying magnetic field. The energy is such 
that a large number of levels are excited, but still low enough such that tunnelling induced 
by an anisotropic perturbation is weak. We shall show that under different conditions (de- 
pending on the fields parameters), different types of magnetization dynamics are realized. 
For the time evolution of the magnetization vector, under certain conditions we obtain a 
solution of the soliton type. The various types of the dynamics will be linked to the structure 
of the phase space of the system. In particular, the existence of the separatrix in the phase 
space has a profound influence on the system behavior. We will show that in this case a 
new type of field-assisted magnetization dynamics emerges, namely a dynamically induced 
switching. This occurs when the energy of the system (in the presence of the field) is still 
lower than the re-scaled anisotropy barrier 1(| and coincides with the separatrix values of 
the energy. Therefore, the domain close to the separatrix is identified as the phase space 
area where the dynamically induced switching takes place. 

II. MODEL 

We consider a molecular magnet, e.g. Fe$ or Mriu acetate. The uniaxial anisotropy axis 
(easy axis) sets the z- direction. The MM is subjected to a constant magnetic field directed 
along the x-axis and a radio frequency (rf ) magnetic field polarized in the x — y-plain. The 



Hamiltonian of the single molecular magnet reads 11] 



H — H + Hi, 

H = -DS2 + gf i B H S x , (1) 
Hi = --gtisHie^iSy + S x ) + H.c. . 

Here D is the longitudinal anisotropy constant, S x , S y , S z are the projections of the spin 
operators along the x,y,z axis, g is the Lande factor, and /zg is the Bohr magneton. H 
stands for the constant magnetic field amplitude whereas Hi, and u>q are the amplitude and 
the frequency of the rf field. The problem when both fields H Q , Hi are time dependent was 



studied in 17] . Using quantum-mechanical perturbation theory, the probability of quantum 
tunnelling of magnetization has been estimated. However, here we are interested in the 
exact solution of the semi-classical equations of motion. Typical values of the parameter D 



are 90 GHz for Mn X2 and D = 30 GHz for Fe^ 
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15| . Since we are interested in the case 



when large number of levels are exited, the spin of the magnetic molecule can be treated as a 
classical vector on the Bloch sphere. Taking into account that S 2 = S 2 +S 2 +S 2 is an integral 
of motion, it is appropriate to switch to the new variables (S z , if) via the transformation 



[14]: S x = a/1 — S z cosy, S y = a/1 — S z siny and rewrite (1) in the compact form: 
A 



H = ~S 2 Z + a/1 - S 2 Z cosy - £a/1 - 5"2(siny + cosy) cos{uj Q t). (2) 

Hereafter, if not otherwise stated the energy and the time scales are set by the constant 
magnetic field H i— > H/gfiBH S, t i— > uj i— > ^s^o- We introduced two dimensionless 

parameters A = g ^^ , £ = ^ < 1- The corresponding Hamilton equations are 

S z = — — — = a/1 — S 2 siny + e\f \ — S 2 (cosy — siny) cos(co>o£); 

<9iJ . cosy . S z . , 

y = — = -(A + = ^ 2 )g z + e = = (smy + cosy) cos(^ t). (3) 

These equations are nonlinear. Therefore, the solutions to (3) can be regular or chaotic, 
depending on the values of the magnetic fields (parameters A, e). From the intuitive point 
of view it is obvious, that for the low energy case, i.e. close to the ground states S z « ±1, 
the system Eq.(2) should become linear. However in the language of variables action angle 
(S z , y) that is not so trivial. Therefore, we will discuss this question in more details when 
studying solutions of the autonomous system. 



III. AUTONOMOUS SYSTEM: AN EXACT SOLUTION 

We inspect at first the autonomous system, i.e. when e = 0. In this case the system can 
be integrated exactly: Taking into energy conservation H = const = — X 

^-A/T^Sfcosy = S (4) 

and 



Sz = a/1 - -S^siny, (5) 

we find 



S*+[^-Z] 2 = l-S 2 . (6) 



Consequently from eq. (6) we infer 



s„(o) 

Xt f dS z 
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s,(t) 



(ir(i-^)-^-fr 



(7) 



This relation can be rewritten in the form 



8.(0) 

dSz (8) 
(a 2 + S 2 ) {b 2 - S 2 ) 
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Sz(t) 



where a 2 = ^ [# 2 /2 - (£A - 1)], b 2 = £[9 2 /2 + (£A - 1)], 9 2 (\) = 2^A 2 -2SA + 1. 
Performing the integration (8) and inverting the result we obtain 

s fbax[(b\/k){t-a), k], < Ar < 1, 

[ bdn[(b\/k)(t-a), 1/k], k>l. 

Here cn(...) and dn(...) are the Jacobi periodic functions. The coefficients that enter eq. (9) 
read 



, 2 If b\ \ 1 



a 



i + < EA - i; 



VA 2 + 1 - 2SA 



2\6(X)J 2 

2[Av / a 2 + 6 2 F(arccos[^(0)/6],A;)] _1 . (10) 



With F((p,k) = J dq(l — k 2 sin 2 l) -1 / 2 being the incomplete elliptical integral of the first 
o 

kind. From eq.(9) we conclude that, depending on the values of the parameter k (10), 
the dynamics of the magnetization is described by different solutions. They are separated 
by the special value k — 1 of the bifurcation parameter k indicating thus the presence 
of topologically distinct solutions. In eq.(9) the Jacobian elliptic functions cn(<p,k) and 
dn(</?, k) are periodic in the argument tp with the period AK(k) and 2K(k) respectively, 



where K(k) = F(ir/2,k) is the complete elliptic integral of the first kind 
period of the oscillation of the magnetization S z (t) is given by 



The time 



*m for < k < 1, 
bX (11) 

2 -W±± for k>\. 



If k — > 1, the period becomes infinite because K{k) — > ln(4/\/l — k 2 ). The evolution in 
this special case is given by the non-oscillatory soliton solution 

S z (t) = b/cosh[b\{t-a)}. (12) 

Considering eq. (10), we infer that the bifurcation value of the parameter k = 1 is connected 
with an initial energy of the system via the ratio 

S 5 = -H s /gfi B H = 1, H(S z (t = 0); <p(t = 0)) = - gi i B H = H s . (13) 



If this condition (13) is not fulfilled the dynamics of the magnetization is described by the 
solutions (9). Finally to conclude this section we consider linear limit of solutions eq.(9): 

cn(u, k) ~ cos(u) + k 2 sin(u)(u — -sin(2u)), k 2 <C 1, 

and 

dn(«, k) « 1 - sin(u) 2 /k 2 , k 2 > 1. 

The interpretation of those asymptotic solutions is clear. First one corresponds to the 
case when in the effective magnetic field H e ff = (g^sifo, 0, —DS Z ), the x- component is 
dominant. Therefore the magnetization vector performs small oscillations ^(t)! < 1 trying 
to be aligned along effective magnetic field. While in the second case, corresponding to the 
ground state solution (system is near to the bottom of double potential well) the effective 
magnetic field is directed along the z- axis. 



IV. TOPOLOGICAL PROPERTIES OF SOLUTIONS 



As established 
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201 ] . the existence of a bifurcation parameter indicates that the solu- 
tions separated by it, have different topological properties. Therefore, it is instructive to 
consider the properties of the solutions (9) in the phase plane. The existence of the inte- 
gral of motion (4) in the autonomous case makes it possible to express S z as a function of 
(f: S z (ip,H). The phase portrait of the system is shown in Fig. ([I]): The different phase 
trajectories correspond to the solutions (9). The phase trajectories corresponding to the 
solution S z (t) = da(ip,k), k > 1 are open and they describe a rotational motion of the 
magnetization. Trajectories corresponding to S z (t) = cn(ip,k), k < 1 are closed and they 
describe the oscillatory motion of the magnetization. Closed and open phase trajectories 
are separated from each other by the special line called separatrix. The existence of a sepa- 
ratrix is insofar important as the states in the phase-space area near the separatrix are very 



sensitive 20| to external perturbations, which signals the onset of chaotic behavior. The 
role of perturbations in our particular case is played by the applied periodic magnetic field. 
We recall that the stochastic layer has finite size and it occupies a small part of phase space. 
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FIG. 1: Two types of phase trajectories of the system separated by the separatrix k = 1, X = Eg. 
The open trajectory (solution S z (t) = dn(i, 1/k), k = 1.52, £ > ) corresponds to the rotational 
regime of motion. The closed trajectory (solution S z (t) = cn(t,k), k = 0.89, £ < ) to the 
oscillatory regime. The separatrix crossing point (^) is of special interest: around this point any 
perturbation leads to the formation of homoclinic structure. 

V. FORMATION OF A STOCHASTIC LAYER 



To determine the width of the stochastic layer we follow Ref. 20|. For details of the 
brmation of the stochastic layer and for the general formalism we refer to the monograph 
2o| . Here we only present the main findings. We introduce the canonical variable of action 
I — § S Z (T>, ip)dip and rewrite the driven nonlinear system (2) in the following form: 

H = H Q + £V(I,<p)cos(u t). (14) 



Here H = u(I)I, u(T) 



-1 

d/(S) ' 



The trajectories laying far from separatrix of the 



unperturbed Hamiltonian Hq are not influenced by perturbation. The motion near the 
homo-clinic points of the separatrix is very slow [2^] . Because the period of motion described 
by (11) is logarithmically divergent, even small perturbations end up with a finite influence 
due to the large period of motion. Thus, the equations of motion for the canonical variables 
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dl ■ e dV ■ , . ,., 

= TTT^TT^COS^ot), (15) 



<9# as* 

<^ = — = +e— S^cos(^), (16) 

may be integrated taking into account the features of the motion near to the separatrix. 
Namely, the acceleration S z gives a nonzero contribution in the integral J dt^S z cos(uj t) 
only near to the homoclinic points 20J (the particle moves along the phase trajectory very 
fast and spends most of the time near the homoclinic points). Therefore, the differential 
equations (15), (16) can be reduced to the following recurrence relations: 

dV ■ 

dt-^-S z cos(u t), (17) 



co(I) J dS z 

At 

P=*>+-7&- (18) 

Here /, Tp, and /, ip are the values of the canonical variables just after and before passing the 
homoclinic point, At is the interval of the time where S z is different from zero. One can 
deduce the coefficient of stochasticity by evaluating the maximal Lyapunov exponent for the 
Jacobian matrix 



dl di_ 

ai d v 

dip dip 

dl dip 



(19) 



of the recurrence relations (17), (18). All of this subsume to the following expression for the 
width of the stochastic layer 

k °-^\jh\- (20) 

Here e, ujq are the amplitude and the frequency of the perturbation. Note that the expression 
(20) is general 2(| and the only thing one has to do is to calculate the nonlinear frequency 
co (I) and its derivative with respect to the energy for the particular system. Thus, even 
for small perturbation (in our case it is the magnetic field with the frequency c^o and the 
amplitude e, see Eq.(2)) the dynamics near the separatrix k = 1, H c = —gfi^HQ is chaotic 
and unpredictable. Consequently, the solutions (9) have no meaning near the separatrix. 
At the same time far from the separatrix H ^ H$, S ^ 1, k ^ 1 they are valid. We note 
that the expression (20) is valid for a low frequency perturbation ujq ^ D and for a high 
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frequency perturbation ujq ^ D as well. For estimation of the width of the stochastic layer 
Kq the variable of action should be determined. Taking into account (4) we find 

1/2 



-^-^ 1 2AE — cos 2 <p ± 2A cos ip J 1 + cos 2 cp 



dip. (21) 



If the static magnetic field is weak then A = g^Ho ^ 1 is a large parameter. In this limit, 
we can simplify expression (21) to obtain 



/(E) = /+(£ > 1) = /-(S > 1) = ^\j^ E [^l) , ( 22 ) 

where E(k) is the complete elliptic integral of the second kind. Taking into account (22) 
the expression for the width of the stochastic layer acquires the following form 

— , K (*) E (*\ (23) 

Condition > 1 of the emergence of stochasticity imposes certain restrictions on the 
parameters of the magnetic field e, u, Ho and on the initial energy S of the system. When 
the energy approaches the separatrix value £ — > 1 the condition K > 1 becomes valid 
even for a very small e <C 1 perturbation. This testifies the fact that the system near the 
separatrix is sensitive to small perturbations. The emergence of chaos is proved by numerical 
calculations as well, see Fig. (J2J). As one can see from this plot, the dynamics is not regular. 
The projection S z (t) of magnetization changes orientation in a chaotic manner. 

However, a chaotic change of orientation is not a reversal to a stationary target state. 
Under dynamical switching we understand here the transition between the oscillatory and 
the rotational types of motion. To be more specific let us discuss the geometrical aspects of 
the motion for the trajectories near the separatrix. Upon applying a static magnetic field, the 
magnetization precessional motion in our case is markedly different from that in the standard 
NMR set up: The key issue is that the effective magnetic field H e ff = (gfj, B H ,0, —DS Z ), 
due to the nonlinearity of the system, depends on the values of S z . The magnetization 
vector tends to align as dictated by the effective field. However, the orientation of effective 
field changes in as much as S z does. Only in the special case S z = 0, ip = 0, 2ir which 
corresponds to the homoclinic points the magnetization vector tends parallel to the effective 



field M\\H e ff. On the other hand, the homoclinic point is an unstable equilibrium point. 
Therefore, the influence of the variable field leads to a switching between the two types of the 
solutions (9). Hence the following scenario emerges: Suppose at the initial time the system 
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FIG. 2: Chaotic motion near the separatrix (k = 1, S = Eg = 1), D = 90GHz. Time independent 
field -ffo) is chosen such that A = — 4, and the ratio between the time independent and 

variable fields is e = H\/Hq = 0.3. The initial energy # = -4.5 • 10 3 GHz is 8/9 of the re-scaled 
barrier height E' B = DS 2 {\ — j) 2 ■ Frequency of the variable field is ujq = 5. One observes that 
the orientation of the magnetization is changing in time chaotically. 

is prepared in the degenerated ground state M s = S. We apply a constat magnetic along 
the x— axis and tune its amplitude to realize the separatrix condition E s = g^H^. A small 
perturbation can then lead to the transitions. In particular, switching off the perturbation 
we end up with the transformed state (cf. Fig. (3)). 




t 



VI. DYNAMICS FAR FROM THE SEPARATRIX: THE MEAN HAMILTONIAN 
METHOD 



To conclude our study, finally we consider dynamics far from the separatrix. The key 
point is the fact that stochasticity emerges in the small phase-space domain located near 
the separatrix. Far from the separatrix the dynamics is regular, even in the presence of 
small perturbations. In this regime, if the frequency of the variable field is high, analytical 
solutions are found with the help of the mean Hamiltonian method. The basic idea of 
the mean Hamiltonian method is the following: for a system having different time scales, 
one averages over the fast variables and obtains thus an explicit expression for the time 
independent averaged Hamiltonian 2l|. In our case, the following condition should then 
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FIG. 3: Motion near the separatrix (k = 1, S = T, s = 1), D = 90GHz, H = -4.5-10 3 GHz, e = 0.3, 
A = 4, ojq = 5. The variable field is applied during the finite time interval between t\ = 100 and 
T2 = 150. Before applying the variable filed, the motion is regular and is of an oscillatory nature. 
The variable field produces a transition into the rotary regime and then is switched off. During 
the transition the motion is chaotic. 



hold: 



u n ( 2DS \ 
gH B H < D < uj {^—) : 



e = Hi/ H < 1. 



(24) 



This condition implies that the amplitude of the magnetic fields should be small and the 
frequency should be high. Provided those conditions hold it is possible to average the dy- 
namic over the fast frequency c^o- The averaged Hamiltonian is determined by the following 
expression: 

- J — 1 1 -~ (25) 



H„ 



H + -{(5H), H} + {(5H), H + l -H}} 



where {A, B} is the Poisson bracket, 6H = H — H, (SH) = J 5Hdt, (. . .) means averaging 
over the time. Applying the procedure (25) to the Hamiltonian (1) and after straightforward 
but laborious calculations with the accuracy up to the second order terms (l/u ) 2 we find 



H av = DS Z + g/j, b H oy /l - S 2 Z cos <p 



{gH B Hi 



x 



X 



^(cos(^) +sm(ip)f 



2^ 2 )(cos(^)-sin(^)) : 



(26) 



The Hamiltonian (26) allows for further simplification: Considering that the variable </? is 
fast in comparison with Si 1 , rotating wave approximation can be used. The Hamiltonian 
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FIG. 4: The dynamics far from the separatrix (k = 1.6, £ = 4Es) is regular; D = 90GHz, 
H = —0.57 • lO^GHz, e = 0.3, A = 100, c^o = 10. The orientation of the magnetization oscillates 
with time, however without a change of sign. Dynamically induced switching is not possible far 
from the separatrix. The left plot corresponds to the numerical solution of eq.(3). The right side 
corresponds to the solution (9) S z (t) = bdn[6A/A;(t — a), 1/k], with re-scaled A constant (27). The 
solutions are in a good agrement with each other. The only difference is the absence of amplitude 
modulation in the analytical approximation. 

obtained in this way is completely identical to (4). This means, that the solutions (9) are 
still valid. The difference is that, the constant A has a different form and depends on the 
parameters of the variable field 



A 



( H igt M B \ 2 \ 2DS 
2u J J gfi b H ' 



(27) 



By comparing the analytical solutions with the results of the numerical integration of the 
system of equations (3) far from the separatrix we verify the validity of our approximations. 

Fig. (4) is for the parameters of the perturbations that are analogous to Fig. (J2J). However, 
unlike Fig.(j2J), where the system is near the separatrix k ~ 1, in the case of Fig. (4) k — 1.6 
which means that the system is far from the separatrix. That is why the dynamics of 
magnetization is periodic in time. The difference, between Fig. 4. and the analytical solution 
(9) is that the amplitude of the oscillations is modulated in time. This observation can be 
explained with the aid of the average Hamiltonian. The point is that the solutions (9) do not 
account for the existence of multiple angles in the average Hamiltonian that were ignored 
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by us. They may lead to the appearance of breathing and amplitude modulations. 



VII. CONCLUSIONS 

We have considered the spin dynamics of a molecular magnet, when the number of the 
involved levels is large. The dynamics of MM driven by variable filed has been studied 

n n 

before [17J. However, in contrast to |17|], the applied fields in our case are quite strong, 
i.e. we are in the strongly nonlinear, non-perturbative regime. The underlying dynamics is 
then treated semi-classically. We showed that the phase space of the system contains two 
domains separated by a separatrix line. The solutions far from the separatrix correspond 
to the rotating and the oscillatory regime, while the separatrix solution is non-oscillating 
and is of a soliton type. The existence of the separatrix is important as the states in the 
domain near to it are extremely sensitive to small perturbations. Therefore, if a variable 
field is applied, instead of a soliton type solutions, the spin dynamics turns chaotic and 
unpredictable. The control parameter is the initial energy of the system. By a proper choice 
of it each type of the dynamic can be realized. The structure of the system's phase space 
is directly related to the possible mechanisms of the magnetization reversal. Namely, if 
the energy is equal to H s = \E B of the re-scaled anisotropy barrier E B = DS 2 (l — 

(the separatrix condition) an external variable field leads to a chaotic change of the 
magnetization orientation. The switching process is random and with the equal probability 
1/2, the system may appear in the new state as well as stay in the old one. The information 
about initial state is lost. This result is different from the case of weak applied fields 
where the dynamics shows a long-term memory of the initial state. 
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